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Abstract 

The two-dimensional 0(3) model at a vacuum angle 9 = n is in- 
\ vestigated. This model has a severe sign problem. By a Wolff cluster 

(Nl ■ algorithm an integer or half-integer topological charge is assigned to 

each cluster. The meron clusters (clusters with half- integer topological 
charge) are used to construct an improved estimator for the correla- 
tion function of two spins at 9 = ix. Only configurations with and 
2 merons contribute to this correlation function. An algorithm, that 
c~| ! generates configurations with only and 2 merons, is constructed and 

numerical simulations at 9 = tt are performed. The numerical results 
indicate the presence of long range correlations at 9 = ir. 



1 Introduction 

The numerical simulation of many interesting physical systems, including in 
particular quantum chromodynamics (QCD) at non-zero chemical potential 
or at non-zero theta vacuum angle, suffers from a complex action problem 
or a sign problem (see and references therein). The numerical simulation 
of these systems is impossible with standard Monte Carlo methods. In this 
article we investigate the two-dimensional 0(3) model at a vacuum angle 
9 = 7i. This model has a severe sign problem. Unlike gauge theories, for 
which no reliable cluster algorithms are available [13 , W\, the 0(3) model 



can be simulated with the help of cluster algorithms and the meron cluster 
concept provides a solution of the sign problem. The results obtained in our 
studies are of interest for the understanding of certain systems in statisti- 
cal mechanics. Haldane e.g. conjectured |10[ that integer and half-integer 
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one-dimensional antiferromagnetic quantum spin chains behave qualitatively 
differently. While integer spin chains have a mass gap, half-integer spin chains 
should be gapless. In the large s limit (2s + 1 being the number of possible 
spin orientations) there is a mapping from the one-dimensional antiferromag- 
netic quantum spin chain on the two-dimensional classical 0(3) model with a 
9 vacuum term . Integer spin chains correspond to 9 = and half- integer 
spin chains correspond to 9 = it. The 0(3) model at 9 = has a mass gap 
in agreement with Haldane's conjecture. On the other hand, Haldane's con- 
jecture together with the (non-rigorous) mapping of spin chains on the 0(3) 
model implies that the mass gap disappears at 9 = it. This corresponds to a 
phase transition governed by the vacuum angle 9. 



The paper is organized as follows. In section [| the Wolff cluster algorithm 
and the sign problem of the two-dimensional 0(3) model at 9 = tt are de- 
scribed. The meron clusters are defined. They are used for the construction 
of an improved estimator for the correlation function of two spins at 9 = n 
(section |^). In order to calculate this correlation function a meron cluster 
algorithm is developed (section f|). The improved estimator together with 
the meron cluster algorithm will be shown to solve the sign problem. The 
correlation function at 9 = tc is calculated numerically. The results of these 
calculations are presented in section |^. 

2 The Wolff cluster algorithm and the 
sign problem of the two-dimensional 0(3) 
model 

The lattice action of the two-dimensional 0(3) model is given by 

S = -(3Y,e x -e y , /3>0. (1) 

The summation is over all pairs (xy) of nearest neighbor lattice sites. The 
three-component vectors of unit length are denoted by e x and called spins. 
The action has a global 0(3) rotation symmetry. 

In numerical simulations of this model one has to generate spin config- 
urations {e} that are distributed with a probability given by the statistical 
weight factor exp(— ^({e})). In the Wolff cluster algorithm || [7j some spins 
are combined and then updated together in one step. For every pair of near- 
est neighbor lattice sites x and y there is a bond variable b xy which only can 
take the values and 1. If b xy = 1, one says that x and y are connected by 
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a bond. The bond variables define the clusters. A cluster is a set of lattice 
sites. Lattice sites connected by bonds with each other belong to the same 
cluster. 

The Wolff cluster algorithm consists of the following steps: 

1. A vector r of unit length is chosen whose direction (Wolff direction) is 
completely random. 

2. To generate the clusters, nearest neighbor sites x and y are connected 
with probability 

p(e x , e y ) = 1 - exp(min[0, -2/3(e x ■ r)(e y ■ r)]). (2) 

(bond activation probability) by a bond. The probability is largest if 
both spins e x and e y are aligned parallel or anti-parallel to the Wolff di- 
rection r. It is zero if the projections of the spins on the Wolff direction 
have different signs. 

3. A cluster is flipped with probability 1/2 independent of all the other 
clusters. Cluster flip means that all the spins belonging to the cluster 
are reflected on the plane perpendicular to the vector r (Wolff plane) 

e x -» e* = e x - 2(e x ■ r)r. (3) 

e* denotes the reflected spin. 

It has been shown that this algorithm fulfills detailed balance and is ergodic 

EE- 

The spins of the two-dimensional 0(3) model define a mapping from the 
lattice with periodic boundary conditions to the sphere S 2 . These mappings 
fall into different classes characterized by an integer valued topological charge 
Q. For this characterization the spins must be interpolated between the 
lattice points. The charge Q counts, how many times the sphere is covered 
by the mapping. The reason for the existence of the integer valued topological 
charge Q is that the homotopy group of continuous mappings from the torus 
T 2 (the lattice) to the sphere S 2 is isomorphic to the group of integers. 

We consider the two-dimensional 0(3) model with an additional topolog- 
ical term in its action 

S -> S - iOQ. (4) 

with 6 real. A complex probability distribution exp(— S) does not make sense. 
For the calculation of expectation values the additional factor exp(i6'Q) must 
be included in the observable 

(F) - (ex P (i0Q)F) e = ((-1) Q F). (5) 
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For the special case 9 = it the additional factor becomes (—1)^. The factor 
(— 1)^ causes a sign problem. In a Monte Carlo calculation the alternating 
sign (—1)^ leads to cancellations between the contributions from different 
configurations to the expectation value. The numerical calculation of the 
expectation value is difficult because it has a large variance if F is one of the 
usual physical observables. 



Using a Wolff cluster algorithm the flip of a cluster changes the topological 
charge from Q to Q*. We can assign to every cluster a cluster charge || 

that is either integer or half-integer. Clusters with half-integer charge are 
called meron clusters. Note that the flip of a meron cluster changes the sign 
(— 1) Q * = (— \y 2q {— l) Q = (— 1)(— 1) Q . This property of the meron clusters 
is exploited for the construction of improved estimators at 6 = n (see next 
section). 



3 An improved estimator for the correlation 
function 

The numerical calculation of the correlation function ((— l)®e x ■ e y ) of two 
spins at 9 = n suffers from large statistical errors because of the sign problem. 
However the meron clusters can be used to construct an improved estimator 
(see eq. (§)) for the correlation function. This construction is based on a 
proposal by U.-J. Wiese. 

Because of the 0(3) rotation symmetry it is sufficient to consider the 
product of the parallel components of the spins (parallel to the Wolff direc- 
tion). There are two possibilities: Either the lattice sites x and y belong to 
the same cluster C or they belong to different clusters 

<(-i) Q e!. e ;> = 

{{-lfele\Q{C x = C y )) + ((-l^ejeje^ + C y )). (7) 

The symbol G( ) is 1 or 0, if the statement in brackets is true or false 
respectively. C x denotes the cluster that contains the site x. The flip of 
a meron cluster changes the sign (—1)^. Therefore when averaging over 
cluster flips the first term only gets a non- vanishing contribution if the cluster 
configuration contains no merons. If x is contained in a meron cluster Ai x 
the flip of the meron cluster changes (—1)^ and e| simultaneously. Therefore 
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when averaging over cluster flips the second term only gets non-canceling 
contributions if x and y are both in meron clusters and there are no other 
meron clusters 



((-i)°44> 



(8) 



((-1) = C y )e(N M = 0)> + ((-l)°e»e BfM, ^ M y )Q(N M = 2)>. 



The factor (— 1)^ can be removed in both terms on the right hand side. 
Every cluster configuration can be thought of as being generated by cluster 
flips from a configuration in which all spins have been on the same side of the 
Wolff plane. For this configuration (— 1) Q = 1, because Q = 0, and e|e| > 0. 
The flip of a non-meron cluster does not change (—1)^ and e|e|. The flip 
of a meron cluster simultaneously changes (— 1)^ and the sign of e|e| in the 
case M. x ^ M. y . Therefore 



With equation (|9|) the sign problem has turned into the problem of finding 
cluster configurations with only or 2 merons. Usually cluster configurations 
contain many merons and configurations with only or 2 merons are quite 
rare. Therefore one aims at constructing an algorithm that generates only 
configurations with and 2 merons. Such an algorithm is described in the 
next section. The improved estimator on the right hand side of equation (||) 
has no sign problem and therefore enables reliable numerical calculations. 

The first term g M only gets non-vanishing contributions from cluster 
configurations in the zero-meron sector while the second term g 2 M only gets 
non-vanishing contributions from the two-meron sector. In the term #om 
the sites x and y must be in the same cluster. Therefore, if the distance 
r = r(x,y) between the lattice sites x and y is much larger than the cluster 
size, there is no contribution to goM( r )- The decay of ^om(^) is known to be 
mainly governed by the cluster size distribution which is exponential |7|] . On 
the other hand in the term (?2M the sites x and y must be in two different 
meron clusters. The two meron clusters can be far apart. Then 5'2m('") 
gets non-vanishing contributions at large distances r (large compared to the 
cluster size). In this way long range correlations can arise. 

In figure [I] the correlation functions at 9 = and at 9 = tt are compared. 
This figure shows the slower decay of correlation function at 9 = tt. 



<(-ir4 e ;> = ' 

(|eU|e(C, = C y )Q{N M = 0)) + (|4 e J|e(M, + M y )Q(N M = 2)) . 



(9) 
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Figure 1: Correlation function at 6 = and at = 7r. L x L = 40 x 40, 
P = 0.1. 

4 A meron cluster algorithm 

Assume that we have a cluster configuration with at most 2 merons. In 
the Wolff cluster algorithm the Wolff direction, used in the definition of the 
clusters, is not fixed. After a change of the Wolff direction, however, all 
bonds become invalid and have to be calculated anew. The new cluster 
configuration likely contains more than 2 merons. 

Thus, a gradual update (that is an update which changes only a few 
bonds in each step) of the allowed cluster configuration is only possible if 
the Wolff direction is fixed. In this case the spins are only reflected on a 
single plane (perpendicular to the Wolff direction). These reflections do not 
change the spin components perpendicular to the Wolff direction. Therefore 
the resulting algorithm is not ergodic. To make it ergodic we must combine 
it with another update step effecting the perpendicular spin components. In 
this additional step only a few spins and bonds should be updated because 
we want to make only small changes on the allowed cluster configuration. 

For the numerical calculations we use an algorithm with two Wolff direc- 
tions. Wolff direction 1 is fixed during the whole simulation. The clusters 
and their charges are defined with respect to this Wolff direction 1. The 



cluster flips on the corresponding Wolff plane are not ergodic. To establish 
ergodicity we use a single cluster algorithm with a Wolff direction 2 that is 
chosen randomly anew for every update step. In one update step a single 
cluster is calculated and flipped with respect to the Wolff direction 2. After 
the single cluster flip the bonds in configuration 1 are calculated anew and 
the meron number constraint for configuration 1 is checked. If the meron 
number constraint is violated the single cluster flip is undone and the old 
bond configuration I is restored (accept reject step). 

After the single cluster flip only those bonds in configuration 1 are cal- 
culated anew whose corresponding spins got changed. Thereby the change 
of the configuration of spins and bonds induced by the update step is re- 
stricted to a relatively small part of the lattice and the acceptance ratio for 
the update step becomes independent of the lattice size. 

With the introduction of a second Wolff direction, the obvious problems 
with ergodicity are removed. There remains however still the posibility that 
independent of the details of the algorithm non-ergodic behaviour could arise 
if certain components of the zero- and two-meron sector are not accessible 
to our algorithm. It could be that configurations in the zero- and two-meron 
sector exist which could be reached only from sectors with more than two 
merons. This would invalidate any approach which deals with selected meron 
sectors only. We have not been able to rule out this possibility analytically. 
However we have not found any dependencies on the initial configuration in 
our results for the correlation function. 



5 Numerical results 

The meron cluster algorithm is used to investigate the large distance behavior 
of the correlation function of two spins at 6 = n. We calculate the projected 
correlation function at 9 = n 




(see figure 0) for various values of the lattice size L = L\ = L2 at fixed 
(3 = 0.1. The spins are summed over one space direction (the 2-direction) in 
order to suppress the contributions from excited states. This corresponds to 
a projection on the zero mode. The data points are fitted with a hyperbolic 
cosine 

h(r) = A cosh ( ( — — r ) ) . (11) 
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Figure 2: Projected correlation function g p and two-meron part of the 
projected correlation function <?2M,p at L = 60 and (3 = 0.1. 



A and m,r are the fit parameters which must be determined. The L de- 
pendence of the masses m n obtained in this way is shown in figure |3|. The 
L dependence of m n seems to be approximately of the form 

C 

m w = —, C = 3.138. (12) 
ij 

The value of the constant C is obtained from a fit. These data suggest 
that in the infinite volume limit the mass m n vanishes and thus long range 
correlations are present. The numerical data confirm the conjecture that the 
correlation function at 9 = tc is long ranged (cf. section §). Also Haldane's 
conjecture about the behavior of quantum spin chains is supported by the 
results (cf. section [l]). 



6 Summary 

An improved estimator has been defined for the correlation function of two 
spins in the two-dimensional 0(3) model at 6 = tt. A meron cluster algorithm 
for the numerical calculation of the improved estimator has been developed. 
The numerical results indicate the occurrence of long range correlations at 
6 = n. 
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Figure 3: Dependence of the mass on the lattice size L at fixed 
[3 = 0.1. The function 3.138/L is also plotted. 
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